/?-Decay Half-Life of the rp-Process Waiting Point Nuclide ^''Mo 
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A half-life of 2.2 ± 0.2 s has been deduced for the ground-state /3 decay of *'*Mo, more than la 
shorter than the previously adopted value. **Mo is an even-even N — Z nucleus lying on the proton 
dripline, created during explosive hydrogen burning in Type I X-ray bursts in the rapid proton 
capture (rp) process. The effect of the measured half-life on rp-process reaction flow is explored. 
Implications on theoretical treatments of nuclear deformation in **''Mo are also discussed. 



I. INTRODUCTION 
A. The Astrophysical rp-Process 

The present model for the rp-process begins with the 
accretion of hydrogen and hehum rich matter onto neu- 
tron stars from nearby companion stars. Gravitational 
energy is released in the form of X-rays as matter reaches 
the intimate regions of the neutron star. The matter 
is compressed as it forms an accretion disk and trav- 
els through the gravitational field gradient towards the 
neutron star, which eventuaUy results in thermonuclear 
burning. The rp-process is a sequence of (p/y) reactions 
and subsequent /3+ decays up to lo'i'.iosrpg J^J. The pro- 
cess begins at ''^Sc, with the seed nuclei produced by 
a series of fast (a,p) and (p,7) reactions on CNO-cycle 
nuclei. The composition of the accreted matter governs 
the overall duration and peak temperature of the nuclear 
burning stage of X-ray bursts [1] . 

The rp-process circumstances prove plausible for the 
synthesis of stable isotopes on the proton-rich side of the 
nuclear chart (p-nuclei) below A = 107. This is because 
the rp-process ashes are finally determined by /3 decay 
from the proton-rich side of /3-stability ■ The observed 
abundances for the p-nuclei ^■^'^'*]VJo and ^^'^^Ru are cur- 
rently at odds with abundance predictions using other 
solar production mechanisms Q . The rp-process can pro- 
duce these nuclei under certain circumstances, sustaining 
X-ray bursts as a possible contributor to their solar abun- 
dances. The amount of rp-process material emitted into 
the interstellar medium depends on the isotopic overpro- 
duction factor, the frequency of burst occurrences that 
produce p-nuclei, and the fraction of material that es- 
capes the gravitational field of the neutron star 4]. It 
has been pointed out that the high gravitational field in 
the vicinity of neutron stars hinders stellar ejection of all 
but the most energetic ash 5] . The bulk of nuclei synthe- 
sized in the rp-process would not escape the gravitational 
field. This limited mass ejection is a major drawback for 



the modeling of X-ray bursts. Nevertheless, the amount 
of ejected material is still being debated and the compo- 
sition of material produced in an rp-process event relies 
on experimental data, and in instances where these data 
are lacking, nuclear structure models 6]. The half-life of 
*'*Mo is a necessary experimental parameter for modeling 
reaction flow above mass 84, and directly determines the 
amount of ^^Sr formed during X-ray bursts. 



B. Shape Deformation Along N = Z 

Deformed nuclei in general exhibit different single- 
particle level spacing and higher level densities than 
spherical nuclei. Single-particle levels, level densities 
and nuclear masses are important ingredients for cal- 
culating proton-capture rates and /3-decay half-lives [7]. 
The ratio of the 4+ and 2+ yrast state energies [-R4/2 = 
E{A^)/ E(2'^)] can be used as an indicator of shape de- 
formation in even-even nuclei, with a smaller value rep- 
resenting a less deformed nucleus The i?4/2 trend 
for even-even nuclei along N = Z reveals a deformation 
maximum at ™Sr and «OZr (see Fig.[I]) A i?4/2 

ratio of 2.52 for ^''Mo marks the beginning of a transition 
towards the believed spherical, doubly magic ^'^"Sn. 

Theoretical predictions for the /3-decay half-life of 
^^Mo within the Quasi-particle Random Phase Approxi- 
mation (QRPA) vary from 2.0 s by Sarrigurcn et al. [13 



to 6.0 s by Biehle et al. [13|. Fig. [5] shows a comparison 
of the half-lives calculated by Sarrigurcn et al. (QRPA- 
Sk3, QRPA-SG2) and Biehle et al. (QRPA-Biehle) with 
experiment. The principle difference between these the- 
oretical treatises is the set of nuclei used to calibrate the 
self-consistent interaction parameters for particle-particle 
coupling strength and nuclear deformation. 

The QRPA-Biehle prediction relies on the nearby nu- 
clei ^*'^°Mo, ^2Ru, and ^"^Pd, which do not exhibit the 
same level (2.1 < R^/^ < 2.23) 0, [3, [H, [H of deforma- 
tion observed in the N = Z region near A = 80 but ap- 
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FIG. 1: Yrast states up to 8^ in even-even N — Z nuclei. The 
transition energies are given in MeV. -R4/2 ratios are shown 
just below 0"'" state. Taken from [lol. ITU. 



proach spherical shapes, lengthening the predicted half- 
lives. The self-consistent parameters for the QRPA-Sk3 
and QRPA-SG2 cases were derived from experimental 
data from nuclei in the region of interest, and mostly re- 
produce the experimental half-lives using self-consistent 
deformations that minimize the energy. The previously 
measured half-life of ^''Mo reported by Kienle ei al. 
falls between the two calculations. This would imply a 
level of deformation unique to the mass region, perhaps 
inconsistent with the observed trend of measured i?4/2 
ratios along even-even N = Z nuclei. 

Recent re-measurements of the half-lives of ^°Zr [l^ 
and 9i.92,93j^j^ shown in Table U Each of these 

measurements are background suppressed by 7-ray gat- 
ing and are s yst ematically lower than those reported by 
Kienle et al. [l8[, at times by more than la. A similarly 
shorter result for a *^Mo re-measurement would validate 
Sarriguren's approach for N = Z nuclei in this region, 
and give consistency with the deformation implied by 
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TABLE I: Experimental half-lives determined by Kienle et at 
[1^ compared with /3-7 correlated measurements [l9l. [20l|. 
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FIG. 2: (color online) Half-lives of even-even N = Z nuclei 
for A = 64 to A = 92 (filled circles) compared to calculations 
using the QRPA. Details of different theoretical self-consistent 
parameters are given in the text The newly measured 

result for ^^Mo is shown as the open circle. 



II. EXPERIMENTAL PROCEDURE 

Nuclei for this study were produced via fragmentation 
of ^^^Xe projectiles accelerated to 140 MeV/nucleon in 
the coupled K500 and K1200 cyclotrons at National Su- 
perconducting Cyclotron Laboratory. The ^^''Xe beam 
was impinged upon a 305 mg/cm^ ^Be target. Frag- 
ments were transported through the A1900 Fragment 
Separator for separation [21]. The A1900 dipole settings 
were Bpi,2 = 2.9493 Tm, and ^^3,4 = 2.5635 Tm with a 
180 mg/cm^ Al wedge placed at the intermediate image. 
Fragments within ± 0.5% of the central momentum were 
transported to the A1900 focal plane. 

The momentum distributions of nuclei produced in 
intermediate-energy fragmentation reactions are asym- 
metric. In Fig. [3] is shown a simulation of yields as a 
function of Bp for the principle isotopes produced in this 
study. The simulation was performed with the program 
LISE ^22*1 . The production rates from the low- momentum 
tails of more stable species exceeds that of the exotic 
^^Mo such that a decay experiment would not be fea- 
sible without additional beam purification beyond that 
achieved with the A1900. 

The new Radio Frequency Fragment Separator (RFFS) 
was implemented at NSCL to purify beams of neutron- 
deficient nuclei [23] . The device consists of two horizon- 
tally parallel plates 1.5 m long, 10 cm wide, and 5 cm 
apart. A time- varying electric field is applied across the 
plates with an amplitude up to 100 kV peak-to-peak. A 
voltage of 47 kV peak-to-peak was applied for this study. 
This field varied sinusoidally at a frequency of 23.145 
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FIG. 3: (color online) Simulation of the yields in particles per 
second (pps) as a function of Bp for the principle isotopes 
created in this study. 



FIG. 4: (color online) Plot of the vertical beam position as a 
function of TOF. The rectangle highlights the region accepted 
through the 10 mm slit setting. 



MHz, matched to the K1200 cyclotron. The time re- 
quired for nuclei in the study to pass through the full 
1.5-m length of the RFFS was roughly equal to 1/3 of a 
RF period. The Time-Of-Flight (TOF) separation of the 
beam species resulted in each fragment seeing a different 
1/3 of the sinusoidal voltage cycle, and therefore receiv- 
ing each a different vertical "kick" . A diagnostic box was 
located 6 m downstream of the RFFS exit and contained 
an adjustable slit system and particle-counting detec- 
tors. Two retractable parallel- plane avalanche counters 
(PPACs) were placed along the beam path, one upstream 
of the slit and the other downstream. The vertical de- 
flection of the fragments at the slit position is presented 
in Fig. 31 The deflection range selected by the 10 mm 
slit gap is also indicated. The phase of the RFFS can be 
delayed relative to the cyclotron RF such that the frag- 
ment of interest appears at a trough, peak, or node of 
the sinusoid as desired. 

Table HIl quantifies the overall and isotope specific re- 
jection factors observed during this study. The table 
demonstrates the importance of distinguishing the de- 
vice's selective rejection rate as opposed to the overall 
rejection rate. Selective rejection of key contaminants 
^*Rb and ^^Kr was of order 10^, while the counting rate 
of the fragment of interest ^''Mo was not affected. The 
overall rejection factor was 180. Fragments that passed 
through the RFFS slit system were delivered to the ex- 
perimental end station for analysis. 

The Beta Counting System (BCS) [24] was used 
for event-by-event correlation of fragment implantations 
with their subsequent /3-decays. The system employed 
a Double-Sided Silicon Strip Detector (DSSD), a single 
995 /im X 4 cm x 4 cm silicon wafer segmented into 40 
1-mm strips in both the x and y dimension, providing 
1600 individual pixels. The high degree of segmentation 
and beam filtering were both necessary to realize an av- 



erage time between implantations into a single pixel that 
exceeded the correlation time, which is set in software 
based on the expected half-life of the isotope of interest. 

The distribution of ^"'Mo fragments on the surface of 
the DSSD is presented in Fig. O As noted above, back- 
ground is controlled in the experiment by taking advan- 
tage of the high degree of segmentation of the DSSD. 
Here the ^''Mo fragments are distributed over ^ 1/4 of 
the DSSD surface, and overlap in position most strongly 
with ions but not ^^Sr ions, where and ^^Sr 
are two of the three significant contaminants in the sec- 
ondary beam. The other sizeable contaminant, ^'^As, will 
have minimal impact on the accurate correlation of ^^Mo 
/3-decay events, since its half-life is of order days (see Ta- 
ble |TT|. The percentage of ^''Mo implantations that were 
preceeded by a ^"Y implantation in the same pixel was 
evaluated for the 10-s correlation time used to establish 
the decay curve for regression analysis as outlined be- 
low. These so-called "back-to-back" implantations oc- 
curred for ~ 8% of the ^''Mo implantation events, and 
were essentially distributed equally in time as the av- 
erage time between implantations in a given pixel well 
exceeded the correlation time. The impact of the "back- 
to-back" ^''Mo-^^Y implantations is further mitigated by 
the fact that ®°Y has two /3-decaying states as listed in 
Table [III The isomeric 1^ state of ^"Y decays with a half- 
life of 4.8 ± 0.3 s, and is a more significant contributor 
to the background analysis of the ^''Mo (3 decay than the 
longer-lived 30.1 ± 0.5 s ground state with J'^ — 4~ [25j . 
(3 decay of the 1~ isomeric state contributed only ~10% 
to the total ®°Y activity, determined from analysis of the 
delayed 7-ray spectrum of ^''Y collected in the present 
work. 

Prompt and delayed 7 rays were monitored at the ex- 
perimental end station with 16 Ge detectors from the 
Segmented Germanium Array (SeGA) [2^. The detec- 



TABLE II; Isotope specific and overall beam rejection 
rates. The rejection factor is reported as the ratio of 
RFFSo///RFFSon for individual fragment yields. Rejection 
factors larger than 1 indicate a reduced transmission with the 
RFFS on. Rates for specific isotopes are given in particles 
per second per particle nanoamperes of beam, normalized to 
the measured rate of **Mo. The beam intensity was attenu- 
ated when RFFS was off (V = kV) to preserve Si detector 
longevity. 
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tors were mounted on a frame designed to closely pack the 
cylindrical Ge crystals in two rings of 8 detectors around 
the BCS chamber for 7-ray detection. Each detector was 
mounted with its cylindrical axis parallel to the beam 
axis. The DSSD was located in the plane that separated 
the upstream and downstream Ge detector rings, thereby 
maximizing the overall detection efficiency of 7-rays emit- 
ted from nuclei implanted in the DSSD. 

Ions implanted in the DSSD were monitored by loca- 
tion for positron emission in the implantation pixel, as 
well as the four nearest-neighbor pixels, over a correla- 
tion time of 10 s. A correlation efficiency for f3 detection 
of ~ 38% was realized. 

Three Si PIN detectors (PINl-3) placed upstream of 
the DSSD served as active degraders. PINl was used for 
the AE and TOF start signal during particle identifica- 
tion (PID) . The active degrader thicknesses were selected 
such that fragments were stopped in the front 1/3 of the 
DSSD, increasing the probability of detecting the small 
AE of a /3-particle emitted in the downstream direction. 

Downstream of the DSSD were six 5 cm x 5 cm Single- 
Sided Silicon Strip Detectors (SSSDl-6), which were pri- 
marily used to veto signals from light particle that were 
transported to the end station. Each SSSD was seg- 
mented into sixteen strips on one face and aligned so that 
the segmentation in each successive detector alternated 
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FIG. 5: (color online) Distribution of ^^Sr, ^''Mo, and *°Y 
fragments in the DSSD. The vertical deflection of fragments 
by the RFFS based on time-of-flight produces a non-uniform 
distribution of fragments in the y dimension. 



in X and y. 

The DSSD provided a trigger for the master gate on 
either an implantation or decay event that registered a 
coincidence signal in the front and back strips. SeGA 
readout was opened for 20 /iS following each event trigger, 
enabling coincidence measurements of fragment-7 and /3- 
7 events. A coincidence register monitored signals from 



5 



the other Si detectors in the BCS telescope, reducing 
dead time by signaUng readout only for energy signals 
within the master gate. 



III. RESULTS AND DISCUSSION 

The decay curve for /3-decay events that occurred 
within 10 s in the same pixel or 4 nearest-neighbor pixels 
of a ^''Mo implantation is shown in Fig. [6l These data 
were fit based on the maximization of a Poisson probabil- 
ity log-likelihood function that considered the exponen- 
tial decay of the ^^Mo parent, the exponential growth and 
decay of the *'*Nb daughter using a T1/2 of 9.5 s [1^, and 
a linear background. The deduced half-life was 2.0 ± 0.4 
s based on a sample size representing 1037 implantations. 
The reliability of the analysis approach is demonstrated 
by the reproduction of the known half-lives of *'^Nb and 
^^Zi. The half-life of ^^Nb was deduced by Kuroyanagi 
et al. [l^l as 4.1 ± 0.3 s. The value deduced from the 
decay curve in Fig. [6] was 3.8 ±0.2 s, consistent with the 
previous result. A half-life of 5.3 ± 0.5 s for *^Zr was ex- 
tracted by Huang et al. [H, [2^ from a fit of a selective 
7p-gated decay curve. The fit of the decay curve for /3 
events correlated with ^^Zr implantations shown in Fig. 
IHresulted in a half- life value of 5.0 ±0.2 s. A decay curve 
for ^^Zr was also generated by gating on known 7 rays in 
the daughter ^^Y with energies 113, 175, and 230 keV, 
that were observed following (5 decay for the first time. A 
fit to this more selective decay curve, taking into account 
the exponential decay of the *^Zr parent and a constant 
background, culminated in a half-life of 4.8 ± 0.5 s, again 
consistent with the the previously reported half-life and 
that deduced here considering only /3-decay events corre- 
lated with ^^Zr implantations. 

One shortcoming of the fits performed with the Pois- 
son probability log-likelihood function is that all decay 
times were considered to be independent, which does 
not accurately describe the probability density of daugh- 
ter decays. Therefore, decay events were also evaluated 
using a Maximum Likelihood analysis (MLH) f30i| that 
considered /3-decay chains of up to 3 members. Analy- 
sis of the ^''Mo decay chain considered events that oc- 
curred within 20 s in the same pixel of a ^^Mo implan- 
tation. Fixed values of 9.5 s and 25.9 min were used for 
the decay time parameters of the daughter and grand- 
daughter, respectively. Background rates specific to each 
individual decay chain were determined by measuring 
the average rate of uncorrelated /3 events on a pixel-by- 
pixel basis. A half-life of 2.2 ± 0.2 s was deduced using 
all 366 extracted chains. This value is consistent with 
the above described Poisson distribution maximum like- 
lihood method, though the approach is inherently better 
because it relies on a probability density function that 
considers the dependencies of the parent, daughter, and 
granddaughter decays, as well as background. 

The new half-life value is more than la below the value 
of 3.7 1q g deduced from 37 implantations by Kienle et 
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FIG. 6: Decay curves for ^''Mo, ^^Nb, and *^Zr. The solid 
and dashed lines are the parent and daughter contributions, 
respectively, to the decay and include a linear background. 
Details of the fitting procedure are given in the text. 
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al. [l^l- The new, shorter value continues the system- 
atic trend observed in Table U relative to the previously 
adopted value. It is also consistent with the above dis- 
cussed deformation assumed in the Sarriguren et al. (l^ 
self-consistent treatment and affirms the trend towards a 
spherical ^"''Sn implied by the measured ^''Mo R4/2 [lol |. 

The initial intent of the half-life remeasurement of 
®^Mo was to gain additional selectivity in the analysis 
of the decay curve by 7-ray gating. The energy spectrum 
for 7 rays found in coincidence with a ^^Mo /3 decay is 
shown in Fig. [71 Unfortunately, no 7 rays were observed 
that could be attributed to the '*'*Nb daughter. Previ- 
ous work has restricted the ground-state spin assignment 
for ^''Nb to 1, 2, or 3 with positive parity. A 1+ ^'^Nb 
ground-state would be consistent with the low 7 intensity 
observed during (3 decay from the ^''Mo 0+ ground-state. 
Additionally, a 1^ assignment would allow direct transi- 
tions into the the 0+ ^''Zr ground state and explain the 
lack of evidence in Fig. [7] for previously observed states 
populated via ^^Nb f3 decay [3l|. Based on the simplis- 
tic assumption that 100% of the *'*Nb /3-decays feed ex- 
cited states in ^*Zr that depopulate through the known 
2^ level at 540.0 keV, we expected to observe a total of 
14(9) counts at this energy in Fig. [T) The absence of this 
540.0 keV transition is consistent with a portion of the 
^^Nb P decay populating the ^^Zr ground state. This fur- 
ther supports a 1+ spin and parity for the ^^Nb ground 
state. 



The implications of the newly measured ®^Mo half-life 
on the rp-process were calculated using a single zone X- 
ray burst model based on ReaclibVl rates provided by 
JINA Reaclib online database [s^]- The abundance with 
respect to burst duration is shown in Fig. [8] for ^^Mo 
(solid lines) and for all A = 84 isobars (dashed line). The 
shaded region covers the range of previously predicted 
half-lives (0.8 s < T1/2 < 6.0 s) given by various models 
[H, [H, [331 . The dot-dashed line represents the yield 
calculated using the experimental upper limit of 4.7 s 
taken from the previously adopted ^'^Mo Ti/2- The order 
of magnitude uncertainty in the final *^Sr abundance is 
reduced to less than a factor 2 with the new half-life. 
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FIG. 7: Spectrum of 7 rays in coincidence with /3-decay events 
occurring in the same and the 4 nearest-neighbor pixels within 
10 s of a ^''Mo implantation. The transition with energy 385.9 
keV represents the decay of the 2^ state in *°Sr fed from ^^Y 
/3-decay, the principle source of background decay events. 



FIG. 8: Impact of ^''Mo half-life on the final A = 84 isobar 
abundance using a single zone X-ray burst calculation. Abun- 
dance is reported as a ratio of (mass fraction) / (mass number). 
Time corresponds to the hydrogen-ignition start time. The 
solid line, bounded above and below with dotted uncertain- 
ties, shows the result using the newly-measured **Mo half life 
(Ti/2 = 2.2±0.2 s). The dashed line corresponds to summed 
abundance of all yl = 84 isobars. The dot-dashed line repre- 
sents the yield calculated using the experimental upper limit 
of 4.7 s taken from the previously adopted **Mo T1/2 . Shaded 
regions highlight the range in abundance based on half-lives 
predicted previously (0.8 s<Ti/2<6.0 s). 

The amount of material processed above ^^Ni can be- 
come significant during X-ray bursts that last longer than 
10 s. In such a scenario, the T1/2 of ®^Mo can have a con- 
siderable impact on the overall mass processing. A previ- 
ous burst simulation Q based on a low a-separation en- 
ergy for ®*Mo provided by the 1992 Finite Range Droplet 
Model predicted a Zr-Nb cycle via the reaction sequence 
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Escape from this cycle is only possible through ^'^Mo f3- 
decay. The newly reported half-life is still longer than 
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the 1.1 s value used for simulations in [7|, leading to a 
more pronounced bottleneck than predicted. The estab- 
lishment of a Zr-Nb cycle would make the impact of the 
*^Mo for some rp-process scenarios considerable in- 
deed. Measuring the a-separation energy of ^'^Mo is crit- 
ical to determining the existence of such a cycle, and the 
real impact that it has on mass processing above A = 84. 

IV. CONCLUSION 

The half-life for the rp-process waiting point nucleus 
^■^Mo has been re-measured to be 2.2 ± 0.2 s, more than 
1(7 shorter than the previously adopted value. This new 
value is in line with the theoretical predictions of Sar- 
riguren et al. of the mid-mass N = Z region consistent 
with a ^''Mo nucleus that begins a shape transition to- 
wards a spherical -'^"'^Sn. ^'^Mo is an important waiting 
point in the rp-process, determining mass abundance at 
and procession above A = 84. A measurement of the 
®'*Mo a-separation energy is critical to determine the full 



impact of ^ Mo on the rp-process mass flow. 7-ray data 
were not sufficient to restrict the ^''Nb ground-state spin 
beyond what was previously reported, though a tentative 
l"*" assignment is favored. 

The experiment discussed in this paper was the first 
to use the NSCL RFFS. The TOF-specific beam rejec- 
tion capability, critical for selectively reducing the overall 
beam rate, is shown to be sufficient to enable fragment-/? 
correlation of neutron-deficient, longer-lived species. 
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